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It is interesting to investigate the OZI-forbidden radiative decays of orthoquarkonium, 
J/tp — > 77T°, 777 and 777' in perturbative QCD. In this work without the approximations 
adopted in literature we carry out a full one-loop calculation which involves integrations of 
4-point and 5-point loop functions. Our numerical results are in agreement with the present 
' data. We also briefly discuss the decays of J/tp — ► 7 + p°, as well as T — > 7r°7, 777, ^'7. 

o 



1 Introduction 



The OZI rulepQ plays an important role in the processes which occur via strong interaction 
_C ■ and in general the concerned calculations are carried out in the framework of perturbative 

QCD. Thus careful studies on such processes where the OZI rule applies can deepen our in- 
sight to the perturbative QCD. On other aspect, we would encounter another serious problem. 
Namely, even though we can accurately calculate the processes at quark-gluon level in terms 
of perturbative QCD, for evaluating the transition matrix elements, the hadronization which is 
related to the non-perturbative QCD, must be dealt with. Because of lack of solid knowledge on 
the non-perturbative QCD, for a whole calculation, one has to employ some model-dependent 
wavefunctions, which may contaminate the theoretical results. Thus to reduce theoretical uncer- 
tainties, the wavefunctions should be well tested and theoretically studied. Then, although the 
results have more or less model-dependence, one can trust that they are comparatively reliable. 
Moreover, as widely discussed, the radiative decays of J/tp may provide an ideal place to look 
for the mysterious glueballs |2j. Since the structure of glueballs is not clear so far, we ignore 
contributions from possible glueball resonances and just concentrate our focus on the direct 
processes which are supposed to correspond to the experimental data. 

Korner et al.[3| investigated the OZI-forbidden radiative decays of orthoquarkonia in pertur- 
bative QCD and their pioneer work was done more than 20 years ago. Since then, the technique 
of calculating loop diagrams has been improved and the knowledge on the wavefunctions of light 
mesons such as n, 77, rf is much enriched. Meanwhile the corresponding experimental mea- 
surements become more precise and make it possible to test our theoretical understanding on 
both the perturbative QCD calculation and the hadron- wavefunctions which are overwhelmingly 
governed by non-perturbative QCD. 

* Supported by National Natural Science Foundation of China. 
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In Korner et al.'s work [3], the authors tactfully dealt with the complicated Feynman in- 
tegrations of four-point and five-point loop functions, i.e. the D- and E- functions [IJ They 
took the weak-binding approximation for both heavy and light mesons, indeed generally this 
approximation is reasonable and much simplifies the calculations. By the approximation, the 
heavy quarks Q, Q (c and c for J/ip) and light quarks q, q in the decay product— light ortho- 
quarkoinum are set to possess equal momenta and be on their mass shells, i.e. pq = pq, p q = pq 
and p 2 Q = m 2 Q , p\ = m\. 

In the weak-binding approximation where q and q have the same momentum and are on mass 
shell, the flavor dependence of the propagator disappears due to the on-shell condition, thus the 
difference of flavors does not manifest in the integration. Since 7r° is of structure -^(dd — uu), r)g 

is of structure -^(uu + dd — 2ss), the contributions from different flavors cancel each other and 

make the total width null, whereas only for t]q = 775 + dd + ss), all contributions are added 

up and the result is non-zero. Therefore to get non-zero widths for J/ip — > 77r°, 77/8, one must 
abandon the weak-binding approximation, namely, consider the relative momentum between q 
and q and not let them be on mass shell. 

Noting the disadvantage of this approximation, Yang ^ [7j dismissed it and re-derived rel- 
evant formulas. In his work where he neglected the mass of the light quark mass of the light 
meson and obtained an improved analytical expression for the rates of J/ip — ► 777, 7/7. 

As J/ip is much heavier than the produced light mesons (tt°, rj, 7/), the 3-momentum of 
the mesons is large, namely much larger than the relative momentum between q and q, the 
equal- momentum approximation as well as the approximate on-shell condition are reasonable. 
However, this approximation probably too simplifies the picture. Concretely, J/ip (or T) is an 
SU(3) singlet, tt° belongs to isospin-1, 77 has a large fraction of 7/8, and in 7/', the 770 component 
dominates. Strictly, the electromagnetic interaction does not demand an isospin conservation, 
i.e. the isospin of photon can be either or 1. But as we take the most probable structure 
that isospin of photon is 0, the J/ip — > 771" is an isospin- violating process and J/ip — ► 777s is an 
SU(3) violating process and only J/ip — > 77/0 is an SU(3) conserving process. One can expect 
the sequence T(J/ip — » 77r°) < T{J/ip — » 77/3) < T{J/ip — > 7770) and the present experimental 
data on J/ip confirm this conjecture [H]. In our calculations, we will show that this sequence 
is completely understood in a full loop calculation, i.e. the cancellation among different quark- 
flavors in 7T° and 773 causes this pattern. 

We are motivated to re-evaluate the processes J/ip — > 771" , 777, 777' by keeping an arbitrary 
relative momentum between q and q and non-degenerate masses of the light quarks (u, d, s). We 
notice that just because of this mass difference the three modes have different non-zero branching 
ratios. In the early work the authors supposed that all processes occur via J/ip — > 770 + 7, 
and the mixing of 77s and 770 results in the non-zero rate for J/ip — > 777. 

However, as dismissing the weak-binding approximation, the advantages for simplifying the 
calculations are lost, namely, one cannot approximately reduce the five-point and four-point 
loop functions into simple three-point loop functions. It is easy to understand that under the 
weak-binding approximation where q and q have the same momentum and are on shell, the 
numerators of the integrands can be properly decomposed into several groups and each of them 
can cancel certain factors of the denominators, so that the number of the Feynman parameters 
can be reduced, whereas without the approximation, generally all the terms of the 5-point and 
4-point loop functions remain and the expressions cannot be further simplified (see below for 
details.). 
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Following the standard procedure £|J, we deal with the 5-point loop functions and then in 
term of the program [5], we evaluate integrations of the 4-point and 3-point loop functions to 
obtain an effective vertex jlUj for J/i/j — > ^ygg — ► 7 + P where P is a pseudoscalar. 

Moreover, there seems to be another contribution from the tree diagram 

cc( J/VO -»■ 7* ->■ 7 + qq(P) 

shown in Fig.l (a) and (b). By a simple analysis, one can immediately note that the contributions 
from (a) and (b) exactly cancel each other as long as the wavefunction of P is symmetric to the 
two light-quark constituents. Therefore in the radiative decays of orthoquarkonia, the leading 
contribution comes from the one-loop OZI forbidden processes jSJ E]- 

To obtain the decay amplitude, one needs to evaluate the hadronic matrix elements 

< p\v eff m >, (i) 

where P stands for 7r°, 77 and rj . In the calculations, we use the light-cone wavefunction 

[HI C3 EU for p. 

In the non-relativistic model, the wavefunction of the light meson P at zero-point Rps(0) is 
responsible for the hadronization effects [31 HUj. 

The same problem has also been investigated by some other authors |14| 115]. especially an 
anomaly for T — ► 777' is addressed. In their work, the gluon contents in 77 and rf are considered 
and non-perturbative effects are taken into account. In the work the anomaly seems to be 
alleviated in this approach, but the concerned assumption is still under suspicion 16 . We will 
come back to this point again in our last section. 

The paper is organized as follows. After this long introduction, in Sec. II, we present our 
formulation and in Sec. Ill, we make the numerical evaluation of the decay rates of J/ip — > 
771" , 777, 77/ and some necessary input parameters are explicitly given. To investigate the whole 
scenario, we further calculate the rates of J ftp — > y0°7 as well as T — > tt°j, 777, 7/7. We find 
an anomaly for T — ► 7/7. The last section is devoted to a simple discussion and our conclusion. 
Some of the tedious details are collected in the appendices. 



2 Formulation of process V — > P7 

2.1 Derivation of the effective operators 

(i) For V -> P7. 

First we derive the formula for J/ip — ► Pj, where P stands as a neutral pseudoscalar meson 
(7T°, rj and 7/). As mentioned above, the tree diagrams shown in Fig.l (a) and (b) make null 
contributions due to a mutual cancellation among them. The corresponding Feynman diagrams 
which are OZI suppressed and offer the leading contributions to the processes, are shown in 
Fig.2 (a) through (f). 

The effective vertex for J/ijj — ► 7 + P where P stands for a pseudoscalar is in the form ]17l 118] 

V eff = g eff e aP , u p a v e^ q ^, (2) 

where Py=Py = Mj/^ is the four- momentum of J/i/j in its center-of-mass frame, and e* u 
are the polarizations of J ftp and the emitted photon, q is the relative momentum of the photon 
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and the pseudoscalar. The g e ff is an effective coupling and should be derived by evaluating 
the corresponding Feynman diagrams and loop integrations. Concretely, this expression is at 
the hadron level, thus our strategy is to derive all the effective operators at the quark level by 
carrying out the loop integrations and then in terms of the wavefunctions of the mesons, we 
obtain the decay amplitude at the hadron level, indeed our goal is to derive g e ff in eq.(2). 

Without the weak-binding approximation which was adopted in literature[3*lllU|. and keeping 
the masses of the light quarks at the propagators, we re-formulate the amplitude. The amplitude 
of V — > P7 can be divided into three pieces which correspond to Fig. 2 (a), (b) and (c) 
respectively. 

For Fig. 2(a), we have 

m\Q Q gjT a T b ®T b T* 
{(\ Pv - q f - ml)(2^ £ [1) 

[Di V (x,m q )(~p v - q) a Ol apv + D p l u {x,m q )m Q Of Lpv + D%(x,m q )(~p v - qfmQ,0* V(7 
+D%(x,m q )m 2 Q Of lu }, 

(3) 



where 



Ofipu = -^9a^gave Paab v Q a p pu Q u q ^ b ^v q + gp^g av £ Paab VQa pl 3UQUq-ib^v q , 

pua 



°%a = -29acrgau£ (3aab VQap fl UQU q <y b 7 5 v q - 2g f3p g au e Paab VQa aa u Q u q 'y b -/ 5 v q 

+9ap,gau£ Paab VQa al 3UQU q -f b ^V q , 

°% = -2igpp:gav^ aah VQ7aUQU<ribl5V q . (4) 

The functions Df \x,m q ), D^x^rriq) are integrals of four-point functions over the internal 
momentum k whose explicit forms are given in appendix A. 
For Fig.2 (b), the amplitude reads 

[$pv - q? - m^](2vr)4 ™ 

[Df (x,m q )(-p v - q) a O) lupv - £>f ' {x,m q )m Q Ol pv - D%(x,m q )(~p v - q) a m Q Ol vrT 
+D2{x,m q )m 2 Q Of lv \, 

(5) 

and the amplitude corresponding to Fig.2 (c) is 



(6) 



where 



0% P a = i£vnpaVQ'ypUQU q 'y l3 "f 5 v q - 2ie vppd VQ^ a UQU q -f e ^v q , 

O^ua = gcag^ Cae VQO- a pU Q U q J e J 5 V q , 

°lpa = gcagp P e Pcae vQ(j a pUQU q "fel5V q . (7) 
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The functions E" pa (x, m q ), E^ix, m q ), E^° (x,m q ) are three independent integrals of five- 
point functions over the internal momentum k and their explicit expressions are also collected 
in appendix A. 

The contributions of the other three diagrams Fig. 2 (d,e,f) are similar to that of the first 
three (a,b,c), to save space, we omit their expressions in the context. 

It is noted that the contribution of Fig. 2 (c) is a five-point Green's function, namely there 
are five propagators in the loop. As long as the weak-binding approximation is adopted, the five- 
point functions E" pa (x,m q ), E^ 1 {x,m q \ E^ a (x,m q ) can be decomposed into sums of two-point 
and three-point functions, so that the calculations are much simplified and analytical expres- 
sions are eventually derived [SUHj- By contrary, without the weak-binding approximation, such 
decomposition is impossible, unfortunately. Following Denner and Dittmaier [3], we decompose 
the five-point functions into a sum of several four-point functions which cannot be integrated 
out analytically. Instead, we will calculate the integrals numerically. Some details about the 
integration are presented in Appendix B. 

(ii) For V -► p°j. 

Serving as a check, we also calculate the branching ratio of J/ip — ► p°7 in comparison with 
that of J /if) — > 7r°7. The Feynman diagrams are the same as that for J/ip — > 7r°7, and the 
explicit expressions of the amplitudes are given in Appendix A for saving space. 

2.2 The hadronic matrix elements 

With the Feynman diagrams of Fig. 2 we derive an effective lagrangian at the quark level. To 
obtain the hadronic matrix elements and then finally the decay rates, one has to evaluate the 
hadronic matrix elements. It is well known that the hadronization happens at the energy scale 
of Aqcd which is in the range of non-perturbative QCD, so far, there is no any reliable way 
to evaluate the hadronic matrix elements. To do the job, we need to invoke concrete models. 
Since the produced meson is relatively light and its three-momentum is larger than Aqcd, the 
light-cone wavefunctions seem to be plausible for description of the light mesons [13 1 1 1 9 | l2fl] . 
The matrix elements are 



where Oj's and Cj's are the operators and their coefficients derived in last sub-section. For the 
pseudoscalars, the SU(3) flavor wavefuctions are 




-ifpPpa / dx(f> p (n,x)^2(Ci(x,m q )Oi)if v E v f3M^ 




uu — dd 



dd + uu + ss 



and 7/8 



dd + uu — 2ss 



V2 



Vo = 



V3 



V6 
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77 and rj are mixutures of 770 and rjs 

r\ = cos 9rjs — sin 6^770, t] = sin 6r)g + cos 9t]q. 

The normalization of the light-cone wavefunction is defined as 

1 

dxcj)(fi, x) = 1. 



The explicit forms of the light-cone wavefunctions of the light mesons can be different. In our 
later calculations, we take three different types of wavefunctions which are given in literatures jl 11 
[pJIUini as 

<t>\{n,x) = 6x(l - x), (8) 
<£ 2 (M,x) = 30x 2 (1-x) 2 , (9) 

3 (^x) = y(l-2x) 2 [l-(l-2x) 2 ]. (10) 
Finally we obtain the hadronic matrix elements of Ma, Mb and Mq as following 

{Pl \M A \V) = iAQ Qgi T°T> » ^^ J5 —^— ? ^ b ) £ dx4ip{li, x) 
{/p/ l 'M,/[-2i<, ( ,„e ewo! >j4 +2if ra „ e p„4] £ (Df(i,m,))(ip v -q)° 

q=u,d,s 

+ifpfv[-^9 a ^9aue l3aab p P b{e v pP v i3 - £ V /3P VP ) 
+g P ^9au^ aab p P b{e va p v (3 - e vP p va )} ( D x'{x,™> q ))<m,Q 

q=u,d,s 

+ifpfv[-^g a a9aue Paab p P b{e v pP v ^ - e v & v p) 

yctPya tycrPya) 

+ga f ,gaue' 3aab Ppb(£yaPyl3 ~ £yf3P v a)} ^ ( D 2( x , m q))(^Py ~ m Q 

q=u,d,s ~ 

+fpfvM v [-2ig^g au e> 3aab p pb e va } £ (D» 2 (x,m q ))m 2 Q }, (11) 

q=u, d,s 



{ p-,\m b \v) = i^tatfri* rt-^, (K _ _ -g ^-( 7) /; ixMv,*) 

{/p/vMv[-2iv«f)^aP?< +2ie^»p w 4] J] (£>r(^.">,))(ij> v -9)° 

q=u,d,s 

-ifpfv{-^g^g a v^ aab p P b{e vP p v p - e v/3 p vP ) 

+gp P gav^ aab p P b{e va Pyp- £ v f3Pya)} ( D l U ( X i m q)) m Q 

q=u,d,s 

-ifpfv[-2gaagaue Paab Ppb(£ v f3P vl M ~ £ vP PvP) 

-2g/3pg a ^ aab p pb (e 
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+ga t igav^ aab P P b(£ v aP v f3 - £ v f3P v a)] Y (-^2 ( x > m <?) ) ( ~~ lY m Q 

q=u,d,s ~ 

+fpfvM v [-2ig^g au e l3aab p pb e va } ]T (D%(x, m q ))m 2 Q }, (12) 

q=u,d,s 

and 

(P<y\M c \V) = m 2 eQ Q gtT a T h ® 3*3*^-^^(7) £ dxfaM 
{fpfvMvlie^papPeyp -2ie u ^ pe p e p s va ] ^ ( E \ P<J { x ^ m q)) 

q=u,d,s 

+ifpfv[gace /3cTe g^p P 9(e vT p vf3 - e vf3 p vT )} ^ ( E 2 (T ( x i m q)) m Q 

q=u,d,s 

+ifpfv[gace (3cTe g fip p P e(e v rP v i3 - e v i3P v r)} Y i E ^i x ^ m g)) m Q}- ( 13 ) 

q=u,d,s 

3 Numerical results 

In this section, we present our numerical results. 

In the numerical computations, there is a mild Infrared (IR) divergence problem. Namely, 
when we carry out the loop integration and convolution integrals of the effective operators with 
the light-cone wavefunction of the produced meson, an IR-divergence emerges, but it is not as 
serious as that in the B-meson decays and can be removed in simple ways. Our strategy to deal 
with the IR problem is standard, namely we assign a small mass to the gluon and vary it to 
check if the result is stable. Practically, we set the small mass to be from 10 -4 MeV to 10~ 6 MeV 
and find that the result has only negligible changes. Therefore we can trust the obtained result 
which is free of the IR problem. Our final results given in all the following tables correspond to 
the gluon mass of 10~ 5 MeV. 

3.1 Input parameters 

The input parameters which we are going to use in the numerical computations are taken 
as follows^ [El El |22j EH]:/j/^ = 551MeV, f r = 7WMeV, f n = lSlMeV, f v = f v ' = 
157MeV, f p = 198MeV, M j/ip = 3096.87MeF, M T = 9460.3MeF, = 134.9766Me^, 
m v = 547.75M e y, m # = 957.78M e y, m p = 775.8MeV, a s {m c ) = 0.26, a s {m b ) = 0.17, 
m c = 1300MeF, mj = 4700Mey, the mixing angle 9 = —11°, and three possible distribution 
amplitudes of pseudoscalar meson are given in Eas. (|8l9ll0|) . 

We will present the the resultant decay rates corresponding to the three different distribution 
amplitudes respectively in the following tables. 

3.2 Numerical results of process V — > P7 

With the above parameters, the theoretical values for the decay width of these processes in the 
rest frame of J/ip, are shown in Table 1. 
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Table 1: The decay branching ratio of J ftp — ► 7T°7, J/ip — > 777 and J/^ — > 77 7 in the rest frame 
of J/t/> and the three columns correspond to the three different parton distribution amplitudes 
of the produced pseudoscalar mesons (tt°, 77, rf). 



Processes 


m u 


777-d 


77l s 


10 5 x BR(<j>i) 


10 5 x BR(</)2) 


10 5 x BR(<p 3 ) 


Experimental data 




1.5 


4 





2.71325 


0.402025 


25.9819 






2 


4 





0.895075 


0.292837 


9.09646 




J/lp — > 7T°7 


3 


5 





1.27478 


0.428238 


11.8012 


(3.9 ±1.3) x 10- 5 




3 


7 





6.43086 


2.21318 


55.5486 






4 


6 





1.64827 


0.565255 


14.1993 






777,„ 


771 ei 


m s 


10 4 x BR(</>!) 


10 4 x BR(<h) 


10 4 x BR((f> 3 ) 






2 


4 


80 


2.41823 


1.68707 


7.23647 






2 


5 


90 


2.46188 


1.69669 


7.44466 




J/V> -► 777 


3 


5 


100 


2.56809 


1.72051 


8.12141 


(8.6 ± 0.8) x 10~ 4 




2 


6 


110 


2.48921 


1.6743 


7.97546 






4 


6 


120 


2.50536 


1.66464 


8.12365 






3 


7 


130 


2.55453 


1.65799 


8.6264 






m u 


m d 


m s 


10 3 x BRifa) 


10 a x BR{(f> 2 ) 


10 3 x BR(<p 3 ) 






2 


4 


80 


1.19297 


1.07802 


1.92765 






2 


5 


90 


1.18068 


1.07066 


1.90548 




J/V -> 77' 7 


3 


5 


100 


1.17033 


1.06154 


1.87308 


(4.3 ±0.3) x 10" 3 




2 


6 


110 


1.15948 


1.05165 


1.84638 






4 


6 


120 


1.14416 


1.03956 


1.8022 






3 


7 


130 


1.12954 


1.02747 


1.76655 





Obviously, the same procedure can be applied to the radiative decays of T —> 7r°(?7, rf)j. We 
calculate the widths with the same parameters by replacing c-quark in J/ip by b-quark in T, 
Qc by Qb and m c by mj,. Then we obtain the numerical results which are shown in Table 2. 
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Table 2: The decay branching ratio of T — > 7r°7, T — > 777 and T — > 77 7 in the rest frame of T 
and the three columns correspond to the three different distribution amplitudes of the produced 
pseudoscalar mesons. 



Processes 


m u 


m d 


m s 


10 5 x RR(<£i) 


10 5 x BR(<h) 


10 5 x BR(<p 3 ) 


Experimental data 




1.5 


4 





0.958428 


0.674172 


2.78981 






2 


4 





0.634743 


0.478474 


1.60618 




T — ► 7T°7 


3 


5 





0.575449 


0.618915 


0.445925 


< 3 x 10~ 5 




3 


7 





1.92777 


2.81325 


0.5723 






4 


6 





0.4889 


0.723972 


0.071549 






m u 


m d 


m s 


10 5 x BRfa) 


10 5 x BR{(f> 2 ) 


10 5 x BR(<p 3 ) 






2 


4 


80 


0.264287 


0.037964 


9.31202 






2 


5 


90 


0.467742 


0.0426104 


13.9632 




T -► 777 


3 


5 


100 


0.772832 


0.0571983 


19.97 


< 2.1 x 10~ 5 




2 


6 


110 


1.20528 


0.0916147 


27.4375 






4 


6 


120 


1.80313 


0.158933 


36.6147 






3 


7 


130 


2.59821 


0.276932 


47.4917 






m u 


"id 


m s 


10 4 x BR{<j> x ) 


10 4 x BR(4> 2 ) 


10 4 x BR( fo) 






2 


4 


80 


9.58491 


8.54894 


14.7284 






2 


5 


90 


9.57026 


8.55094 


14.6 




T -► 7/7 


3 


5 


100 


9.54632 


8.5466 


14.4558 


< 1.6 x 10~ 5 




2 


6 


110 


9.51302 


7.99717 


14.2925 






4 


6 


120 


9.47596 


8.51945 


14.1333 






3 


7 


130 


9.42768 


8.49574 


13.9498 





As discussed above, as a check, we evaluate the decay width of J/ip — > p°j which corresponds 
to an effective three-vector vertex and has not been observed yet. The results are tabulated in 
Table 3. 

Table 3: The decay branching ratio of J/ip — > p°7 in the rest frame of J /if), and the three 
columns correspond to the three different distribution amplitudes of p°. 



Processes 


m u 


m d 


m s 


10 13 x BR{(f)i) 


10 13 x BR(<p 2 ) 


10 13 x BR(<p 3 ) 


Experimental data 


J/V - p u 7 


4 


6 





4.74477 


2.07037 


37.9532 
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4 Conclusion and Discussion 



In this work, we re-study the OZI forbidden radiative decays of othoquarkonia in the framework 
of perturbative QCD. In the process, we do not take the weak-binding approximation or set 
the light quark mass to be zero and carry out a complete integration of the five- and four-point 
functions. In this scenario, we can take into account the SU(3) and isospin violation which would 
result in non-zero rates for J/ip — > tt°j and J/ip — > r/s7- However, on the other side, without 
such approximation, we cannot derive elegant analytical expressions for the amplitudes as done 
by Korner et al. [3] and Yang [Bj. Instead, we need to invoke complicated computer programs 
and we follow the recipe given by Denner and Dittmaier 0] to reduce the five-point functions 
into a sum of four-point functions and then use the developed computer program "LoopTools" 
to carry out the integrations. 

Another difficult point is to evaluate the hadronic matrix elements which are fully governed 
by the non-perturbative QCD. The initial heavy quarkonium (J/ip or T) are composed of only 
heavy quarks, so that the on-shell approximation is reasonable and its contribution to the 
amplitude can be described by its wavefunction at origin, i.e. i?y(0) [HIE]- As we abandon the 
weak-binding approximation for the produced light meson, the simple description in [3] for the 
produced meson by i?P5(0)|31E] can no longer be adopted. Since the produced pseudoscalar 
meson is light, the light-cone wavefunctions seem to be applicable for the calculations. There 
are several typical different light-cone distribution amplitudes for the light mesons, so far, one 
cannot determine which one is the most suitable. Thus we adopt all the three for our calculations 
and the corresponding results are listed in the tables of last section. 

Our results for J/ip — > P7 where P stands for 7r°, rj, rj' , are quantitatively in agreement with 
the experimental measurements. For J/ip — > 7r°7 which is an isospin violating process because 
it is proportional to an asymmetry of u and d quarks, within reasonable ranges of the masses 
of u and d quarks, all the three distribution amplitudes can result in values in agreement with 
data. For J/ip — > 777, it seems that only </>3 can give the values in good agreement with data. 
For J/ip — > 77' '7, all the obtained values are slightly smaller than the data. This small declination 
can be understood, because theoretically the results depend on the masses of light quarks and 
the QCD coupling, (here we do not include the running of a s as in 3J, and experimentally, 
all the concerned widths are small and certain measurement errors are unavoidable. We can 
expect that the CLEO and BES III which will be running in 2007, can make more precise 
measurements to testify the results. Moreover, we calculate the decay width of J/tp — ► p°7 
which has not been measured yet. This result will be measured in near future and the data can 
provide some information about our understanding of both perturbative and non-perturbative 
QCD (the ansatz for light-cone wavefunctions of mesons). 

As we turn to the radiative decays of T, the situation seems peculiar. It is noted that 
T — ► 7r°7 is related to isospin violation, T — > 7/7 is related to SU(3) violation because 77 has 
a large fraction of 77s, only T — * 7/7 conserves SU(3), as rf contains mainly 7/0. Thus one can 
expect r(T — ► 7r°7) < T(T — > 7/7) < T(T — ► 7/7). For radiative decays of J/ip, this sequence 
obviously holds and our calculations confirm this pattern. However, for T, the measurements 
seem not to follow the pattern [H]. Of course the data only set upper bounds on these decay 
modes, there is still some possibility to upset this pattern. 

Our theoretical results for T still follow the sequence, and for T — > 7r°7 and T — ► 777, the 
values are consistent with the experimental upper bonds, but for T — > 77' 7, the calculated value 
is obviously larger than the upper bound set by the present experimental measurement. It may 
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indicate some anomaly in these decay modes of T and the reason is worth further studies both 
theoretically and experimentally. 

The anomaly problem has been observed by several authors [HI 1141 1151 116j , all agree that 
in the perturbative framework which is the main content of this work, this anomaly exists. To 
reconcile the theoretical results with the data, the gluon contents in rj and rj' are considered 
[1411151 HE] , and the non-perturbative matrix elements < 0\GG\rj( > > was evaluated by Feldman 
and Kroll [Uj. It seems that the non-perturbative effects alleviate the discrepancy between the 
theoretical results and data. However, as pointed in ^H], the problem is not fully solved yet and 
needs to be investigated further. In our work, we are only dealing with the perturbative part 
and confirm existence of the anomaly as discovered before. We will investigate the anomaly and 
the probable non-perturbative effects in a wider range. Moreover, in the works [6} I14 | I15 [ IT6] . the 
77T° case was dismissed because of its smallness, by contraries, we explicitly keep the light-quark 
masses and obtain the corresponding results for 771" final state in the radiative decays of both 
J ftp and T. 

Our starting point is indeed the same as that of Korner et al. [H] and Yang pj, except we 
do not take any approximation and carry out the full integration of the five- and four-point 
functions. We find that the numerical results for J /ip are consistent with data, but there is an 
anomaly for T —* 7/7 if the present measurement is correct. All these need more and careful 
investigations. 
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Appendix A 

The integrations corresponding to Fig. 2 (a), (b) and (c) are 
Ma = I -^e^(i)v Q {p 2 )(-ig s T a la )- — 1- {-ig s T b lp )- 



(27r) 4 j>\ - (fi + # - m Q i>\ — i~ ttiq 

(-ieQQj l t)uQ(pi)u q (p3)(-ig s T b ~fP) % (-ig s T a -f a )v q (p A ) 

p3 1 r m q 



—1 



(pi + P2 + k - q) 2 k 2 

= J -^eQ Q g\T a T h ® T^TVgfe) 7^1 - 4 + j& + mgh^i - 4 + m Q ) 7lM u Q ( Pl ) 

1 



k 2 (pi +P2 + k - q) 2 [(pi - q + k) 2 - m 2 ][(pi - q) 2 - m 2 }[{p^ + k) 2 - m 2 ] 



Mb =J ^^(7)^(p 2 )(-ieQ Q 7,) — L (-ig s T" la ) 



(Al), 



(27r) 4 <A - fa - m Q fa+ft-mQ 
{-ig s T b lp )u Q (jp x )u q ( Pz ){-ig s T b ^) 1 —l (-ig s T a 7 a )v q (p 4 ) 

P3 1 ft ~ m q 
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(pi + Pi + k — q) 2 k 2 

= J ^^(7)eQ c ^T a T 6 ®r fe rVQ(p 2 ) 7 ^(^-^ + mQ) 7Q ( l 5 1 + ^ + mQ) 7 ^Q(p 1 ) 

Uqipsh^ifa + ft + m q )j a V q (p 4 ) 

\ ^2) 

k 2 (p!+p 2 + k- q) 2 [(q -p 2 ) 2 - m 2 Q ][(p 1 + k) 2 - m 2 Q ][(p 3 + k) 2 - m 2 }' 

d 4 k 



/a Ik i 
_^ (7) ,- Q(p2)( - ifo r» 7 ,) - _ i+ (-feOar,) 



-ig s T h -fp)u Q {pi) 



fii + ft-niQ 



u q (p 3 )(-ig s TV) , , ' (-i9sT a r)v q (p4) 

P3 + ft — 



—I —I 



(Pi + P2 + k - q) 2 k 2 

d 4 k 

(2tt)' 

Uq(p 3 h P (h + ft + m q h a v q ( P 4)e*»W 

(AS). 



= eQ Q g A s T a T b T b T a J JA*r Q ( P2 ) 7Q (^ _ ^ + ft + mQ ) 7/i (^ + £ + m Q ) 7 „UQ( Pl ) 



/c 2 (pi +p 2 + k - q) 2 [(pi -q + k) 2 - m 2 Q \[(pi + k) 2 - m^][(p 3 + /c) 2 - m 2 ] ' 
The expressions for the process J/tp — > p° 7 are 



Z7T 



2 eQ Q g 4 T a T b ®T b T a 



Ma = f(i % 2^I £ * M(7) 



+D»(x,mq)m Q Oi llu }, 



where 



°lapv = ^guaVQ-f^UQUq-fpVq + 2g ap VQ-f u U Q U q -f p V q - 2g pv V Q ^ a U Q Uq^ p V q 
+ C 2g p vVQ^ ll UQUq~f (T Vq + 2g pv g ap VQ~f f5 UQU q ^ V q - 2g pu VQ<y a UQUq<y p V q 
O^pu = 2iVQG vp UQU ( fi rho Vq + 2ig pv VQG fipUQUq-i 13 Vq 
Oflw = 2iVQO ap UQUq~i v Vq 
0% = ZVQ-f^UQUq^Vq. 



((K-?) 2 -^)(2^ ™ 
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[-D^(x,m q )(-p v - q) a Ol apu + D^{x,m q )m Q Of ipv + D%(x, m q )(-p v - q) a m Q 0^ ua 
-D v 2 (x,m q )m 2 Q 0* v }. 

The corresponding hadronic matrix elements are 

(v'\qj a qJ2(Ci(x,rn q )Ol)Q^Q\vj 

= rn , [ dxcf) ,(fj,,x)y2(Ci(x,m q )Oi)if v e v pM v 

V v a V Jq V — 

(v)q la qJ2(Ci(x,m q )Ol)Qap p Q\V^ 

= i K lE v'J n v l f dx ^ v '^ x )^^ C ^ x ^ m Q)°i) i fv i ( £ vPPvp-PvP £ vp)- 

Jo i 

Concretely the hadronic matrix elements of Ma and Mb are 
(V'j\M A \V} 

= i7r 2 eQQggT a T b * 2 ^^ ?J _^__^( 7 ) £ d x M », x) 
{if v , m v , if v M v i 2 9ua£ v ^ v , p + 2g api £ vV e v , 2g iXV £ ytT £ ^ + ' 2 9pv e vt J - e v i a 
+2g pu g afl e vl3 e /3 v/ - 2g pv e v<T e v ,J ^ ( D T ( x , m q))(\p v ~ 

q=u,d,s 

+if v ,m v ,if v i[2i(£ vV p vlx - p vV £ vp )£^ ^ + 2ig pu (e vf3 p vfl - p v p£ vfl )e^} ^ (£>f '(x,m q ))m Q 

q=u,d,s 

+if v/ m v ,if v i[2i(e vr7 p vfl -p va e vfl )e v/ J ( D 2 (x, m q ))(^p v - q) a m Q 

q=u,d,s 

+if v ,m v ,if v M v [2e vP e v ,} ( D 2{x,m q ))m 2 Q ,} 

q=u,d,s 

and 

(V-i\M B \V) 

{-if. m , i/ v M v [2gwe vAt e , 

V V V p V p V p V <7 

+2g pu g a ^e vl 3£ l3 v , -2g pu e va e v ,J J2 ( D i"( x > m q))i\p v ~ <lf 

q=u,d,s 

+if v ,m v ,if v i[2i(£ vl/ p vpj - p v y£ vP )£ v , p + 2ig pu {£ vf5 p vp - p v fs£ vP )£^,] W ( x i m q)) m Q 

q=u,d,s 

+if vl m v ,if v i[2i(e v(T p vp -p v(T £ vP )£ vl } ^ (D%{x,m q ))(-p v - q) a m Q 

q=u,d,s 
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if v ,m ,if v M v [2e VIJ ,e ] ( D 2(x,m q ))m 2 Q .} 



1 ' V' V' " v v L v f~ V' 

q=u,d,: 



Appendix B 

The two four-point loop functions in our text are: 



i r+, (p 1 + k)p(p 2 + ky 



D l ( X > m *) in 2 J d \2 [ ( pi+k) 2_ m 2 Q]{k+p2 y [ ( k+p3) 2_ m 2y 

m l) = -2j d " k k 2 [{pi + k) 2 _ m 2p + l 2) 2 [{k + ps) 2 _ m 2] ' ( 14 ) 



with 



Pi = ^P V ,P2 =P V ~q,P3 =P P x. (15) 
And the three five-point loop functions are: 



in 2 J k 2 [(k + P3) 2 — m 2 ](k + P2) 2 [{k + pi) 2 — rriQ][(k + pi) 2 — m 2 1 ' 



E^(x,m q ) = 1- ( d 4 k 



{k+Pi) u {k + p z y 

k 2 [(k +p 3 ) 2 - m 2 }(k + p 2 ) 2 [{k + P4) 2 - m 2 Q ][(k +p\) 2 - m 2 Q ] ' 



,p°^^x- 1 t {k+Piy{k + p 3 y 



E 3 (X,m q ) .^Jd ^2[( fc + k) 2 _ m 2] (jt + ?2)2 p + p4 )2 _ m 2 ][(Hpi) 2 _ m 2] ( 16 ) 

with 

Pi = 2 p v p2 =Pv -1,P3= P p x,Pa = -p v - q, (17) 
which can be decomposed into somes of four-point loop functions according to 4 . 

+(g ||»"*(.) + El-^W-BK 
+ (g^« </, " , "(^t(-')' +j ^^„^ (i ,) pW 

14 



with 



2 — U — 1 



(y)ij = m- + m) - (p- - Pj) 2 , i, j = 0, . . . , 4, 

(i (4) )« = 2j4p{, fe,Z = l,...,4, 

m = 0, mi = m q , m 2 = 0,m 3 = m Q ,m 4 = uiq, 

Po = °> Pi = P3, P2 = P4, P3 = P2,P4 = Pi- ( 19 ) 

And D^ in ^ vp {iys denote the ultraviolet-finite four-point functions that are obtained by removing 
the ith propagator in the five-point functions. Let us take D^ in ^ vp {l) as an example, and it is: 

D (fin)up^ = _L f d 4 k (20) 

m 2 } k 2 (k +p 2 ) 2 [(k + J> 4 ) 2 - m 2 Q \[(k + pi) 2 - m 2 Q \ 

where only four factors exist at the denominator. 

Yi is obtained from the 5-dimensional Cayley matrix Y by replacing all entries in the ith column 
with 1 and the 3-dimensional matrices Z^ result from the 4-dimensional Gram matrix by 
discarding the ith row and jth column. 
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